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Sliding Mode Control Continued 


We generalize the sliding mode control examples of the last lecture to 
the class of systems 


9 = fa(y,w) 


dics Ape yae COs, “ 


where w € IR’, u € RP, € R"?. The uncertain terms are 6(y,w) 
and the p x p matrix G(7,w), assumed to be diagonal with entries 


Let $(7) be a virtual control law for w that stabilizes the origin of the 
y-subsystem, 7 = fa(y,¢(7)). To drive the trajectories to the sliding 
surface w = (7), we note that s := w — (7) € RP satisfies 


$= fol.) — ED fy (p00) + olyn00) + G(r) 
and let 
w= 6% y,.0) | fol) — “ED fun c0)| +0 


where G(y,w) is a nominal model for G(7,w), and v is to be de- 
signed. Then, 


§= (1-66) |fil.c2) - SP fu(y,00)] + 50,0) +6( qo), 


———— 


=: A(7,w) 


~~ 


which means that the ith entry of s satisfies 
§, =Ai(y,w)+ei(y,w)v;, i=1,---,p. 
We let 
vj = —(pily,w) + po)sgn(si), po > 0, (2) 


where 0;(7,w) is a function such that 


Ai(nw)l 


gi(y, w) pi(1], a). 
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Then the Lyapunov function Vj = is? satisfies V; < —V/2p0 goV/ e 
which guarantees finite time convergence of s; to 0, as discussed in 
Lecture 22. 


Thus the trajectories reach the sliding surface w = #(7) in finite time 
and, if the subsystem 7 = fa(1,¢(7) +s) is ISS with respect to s, then 
4 remains bounded during the reaching phase and converges to zero 
asymptotically during the sliding phase. 


Continuous Approximation of Sliding Mode Control 
To avoid the chattering phenomenon discussed in the previous lec- 
ture, we can employ the continuous function 


ore x/e when x € [—¢, €] 
. 7 sgn(x) otherwise, 


which approximates sgn(-) when ¢ > 0 is a small constant. 


If we implement (2) above with o;(s;) instead of sgn(s;), the Lya- 
punov analysis is unchanged when |s;| > ¢, where the two functions 
are identical. Thus, |s;| > ¢ implies Vj < —V/2p0g0V/ 2 < 0, from 
which we conclude that s; reaches the interval [—e, ¢] in finite time 


and remains in it thereafter. Likewise, if 7 = fa(y,@(7) +s) is ISS 
with respect to s, then 7 converges to a residual set around 4 = 0 
whose size shrinks as e + 0. 


Therefore, the continuous approximation eliminates chattering, but 
guarantees convergence to a small set around the origin rather than 
to the origin. 


Example: For the system 
Xy = X1XQ 


Oxt+u, |0| <2, 


x2 

the virtual control $(x1) = —xf and the variable s := x2 — $(x1) = 
Xo + s result in 

xy = =e + X48, 
which is ISS with respect to s. To drive s to zero we note that 

$= Qxixs + Ox, + 

and apply the control 

w= 2x as ko, 
which guarantees global asymptotic stability of the origin (x1,x2) = 
(0,0) with the discontinuous feedback v = —(2x+ + p0)sgn(s). 


If we apply the continuous approximation v = —(2x7 + p9)0¢(s) we 
achieve convergence to a set which shrinks to the origin as e — 0. 
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Tracking Control 


Consider a model represented in the normal form for input-output 


linearization: 
Z2 = folz,¢) 
f& = fs 
a = Cr 
Cr = b(z,f) +a(z,f)u 
y = C1 


where a(z,¢) and b(z,¢) are imperfectly known, but 
a(z,C) > go >0 


with some positive constant go. In addition we assume the zero dy- 
namics subsystem z = fo(z,¢) is ISS with respect to C. 


This system is of the general form (1) with 7 = [z’,@1,--+ ,f,—1]? and 
w = ¢;, and we can design a virtual control 


Or = —kp_10r-1 — +++ — kt (3) 


with coefficients k,,,---k, such that 


. : ‘ ; : 
= (4) 
1 
Cr-1 hi =k) [Grn 
= Ao 


is asymptotically stable. 


The dynamics restricted to the sliding surface (3) consist of the sub- 
system (4) driving the ISS zero dynamics; therefore the trajectories 
converge to the origin. Finite time convergence to the surface is 
achieved with the standard design approach discussed on page 1. 


When the goal is to ensure that the output ¢; tracks the desired tra- 
jectory yq(t), we define the tracking error variables 


e1:= C1 —yalt), 2 = C2 —galt), + er = Gr — y(t) 
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and rewrite the system equations as 


folzer +yalt),..-,¢r-+y¥ (2)) 


Zz — 
& = & 
ert = @, 


é = b(z,0)— y(t) +a(z, C)u. 
As the sliding surface we select 
S:= ey +kp_yep_y +++: +key = 0, (5) 


where k,_1,--++k, are such that the matrix Ao defined in (4) has all 
eigenvalues with negative real parts. Thus, e(t) — 0 on the sliding 
surface and z(t) remains bounded by the ISS assumption, and by the 
boundedness of y_(t) and its derivatives. 


For the reaching phase we note that 


§ = 6 +ky1ép1 +++: +kyé4 
= b(z,f)- y(t) t+ kp_yep +++ thyeo +.a(z,l)u, 
and select 
H= a [EO —IPO+kart the +e © 

This yields 

§= (1- FEB) [1+ G0) - 8.0) +062 00 

a(z,¢) 
; =: A(z,¢,t) 

where [ --- ] is the square bracketed term in (6), and A(z,€,t) de- 


pends on t due to the derivatives of yg(f) occuring in this expression. 


We then choose p(z,Z,t) such that 


|A(z, 6, £)| 
alz,f) < p(z,¢,t) 
and complete the design (6) with 
v = —(0(2,6,t) + po)sgn(s), po > 0. (7) 


Note that, if we set y(t) = 0, the tracking controller (6)-(7) reduces to 
a stabilizing controller for the origin (z,¢) = 0. 


Example: Consider the system 
xX, = X2+sin x1 
Xo = Ox +(140.)u |0,| <2, || < 0.5, (8) 


Y = Xj. 
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To design a tracking controller we first bring the system to the nor- 
mal form with the new variables €, = x; and @7 = x2 + sin x): 


- (9) 
Co = (xo + sin x1) cos xy + 6,x1 + (14 65)u. 


Then we define the error variables e; = €, — yq(t) and e2 = Co — yq(t), 
which are governed by 


ey = 


éy = (x2 +sinx,) cos x, — q(t) + x7 + (1+ 62)u, 
and select the sliding surface 
S:=eo7 +key =0, ky >0. 


Thus, 


§ = (x2 + sin x,) cos x1 — tjqg(t) +kye2 4 0,x7 + (1+ 62)u 
and the feedback 


uo = —(x2+sinx,)cosx, + q(t) —kje2+v 
—(0(x,#) +0) po >9 


S) 
I 


results in 


$ = 02(—(x2 + sin x1) cos x1 + a(t) — kyer) + 4x4 +(1 + 62)v. 
—— 
=: A(x1, X2,1) 


Using the bounds |6;| < 2, |@2| < 0.5 we get 


|A(x1,X2,t)| 0.5|(x2 + sin x1) cos x1 — q(t) + kyen| + 2x7 
1+ 62 — 0.5 
= |(x2 + sin x1) cos x1 Hat) } kye>| } 4x? 


and, substituting e2 = C2 — yq(t) = x2 + sin x; — q(t), we select 


p(x,t) = |(x2 + sinx1) cos x1 — fa(t) + ki (x2 + sinxy — ya(t)| + 4x7. 


It is important to note that sliding mode control can address only 
limited forms of uncertainty. In the example (8) the uncertain terms 
appear in the same equation as the control input; that is, they are 
“matched" to the input. The first equation in (8) contains no uncer- 
tainty, which allowed us to bring the system to the normal form (9). 
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